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ABSTRACT 

We report evidence of a fully established galaxy cluster at z = 2.07, consisting of a ~ 20cr overdensity of red, compact spheroidal 
galaxies spatially coinciding with extended X-ray emission detected with XMM-Newton. We use VLT VIMOS and FORS2 spectra 
and deep Subaru, VLT and Spitzer imaging to estimate the redshift of the structure from a prominent z = 2.07 spectroscopic redshift 
spike of emission-line galaxies, concordant with the accurate 12-band photometric redshifts of the red galaxies. Using NICMOS and 
Keck AO observations, we find that the red galaxies have elliptical morphologies and compact cores. While they do not form a tight 
red sequence, their colours are consistent with that of a > 1.3 Gyr population observed at z ~ 2.1. From an X-ray luminosity of 
7.2 x 10 43 erg s _1 and the stellar mass content of the red galaxy population, we estimate a halo mass of 5.3-8xl0 13 M o , comparable 
to the nearby Virgo cluster. These properties imply that this structure could be the most distant, mature cluster known to date and that 
X-ray luminous, elliptical-dominated clusters are already forming at substantially earlier epochs than previously known. 

Key words. Galaxies :clusters: general - Galaxies: clusters: individual :CL J1449-0856 - Galaxies :high-redshift - large-scale structure 
of the Universe 



1. Introduction 

Massive clusters are rare structures in the distant Universe, 
arising from the gravitational collapse of the highest den- 
sity peaks in the primordial spectrum of density fluctuations 
dPeebles 19931 IColes & Lucchin 19951 IPeacock 1999b . Their 
abundance reflects the original state of the matter density 
field and depends on fundamental cosmological parameters, 
such as the sh ape and normalisation of the matter p ower 
spectrum (e.g. IPress & S chechte rT9741 [Haiman et al. 2001 , 
ISchuecker et al. 2 003 ). Measuring the distribution of galaxy 
clusters can thus place constraints on these cosmological 
parameters and p rovide a powerful tes t o f primordial non- 
Gaussianities JJimenez & Verde 20091 |Cay6n et al. 2010 
ISchuecker et al. 2002 ). As the largest and most massive 
bound structures in the Universe, galaxy clusters are also 
the most biased environment for galaxy evolution and 
constitute a prime laboratory for studying the physical 
processes resp onsible for the formation and evolution of 
galaxies (e.g. IBoselli & Gavazzi 20061 |Park & Hwang 2009 



Demarc o et al. 2010b . The strong dependence of galaxy activity 
on the surrounding environment, which gives rise to the well- 
known correlations of morphological type JPostman et al. 2 005 , 
|Hwang & Park 20091 Ivan der Wei et al. 2010b and decreas- 



ing star formatio n (Hashimoto et al. 1998, IPatel et al. 2 009, 
IRettura et al. 20101 Peng et al. 2010| ) with increasing galaxy 
density, is most easily and dramatically illustrated in the 
extremely dense cores of local massive galaxy clusters. 

But while, about 13.7 Gyr after the Big Bang, today's 
galaxies, baryons and dark matter continue to steadily fall into 
the massive clusters' potential wells, the elliptical galaxies 
that dominate their cores ceased forming stars early on and 
have been evolving passively for most of cosmological his- 
tory. The traces of the formation process of present cluster 
ellipticals thus smeared out, their co-evolution with the cluster 
cores and the assembly history of the latter can hardly be 
reconstructed from low redshift data alone. Furthermore, the 
thermodynamical properties of the baryons in the intergalactic 
medium of clusters (ICM) suggest that ~ 1 keV more energy 
per baryon was injected into the ICM than can be accounted for 
by pure gravitational collapse (Ponma n et al. 1999b . This excess 
entropy was generated presumably by galactic winds powered 
by either supernovae or AGN, but the exact processes have 
not yet been identified. It is also unclear when the correlation 
between the mass, X-ray luminosity and temperature of the 
ICM, crucial for the use of clusters as cosmological tools, were 
first in place. The question of the assembly of clusters, the 
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settling and thermodynamical evolution of the X-ray shining 
ICM within their deep potential wells and the build-up of 
their constituent galaxy population must then be addressed by 
looking as closely as possible at the early stages o f their for- 
mation (IVoit 20051 iPonman et al. 19991 iRosati et al. 20021 ). This 
formative epoch is often put at z > 2, as suppo rted by evidence 
foun d in rece nt years of increased act i vity (Elbaz et al. 20071 
Hayashi et al. 2010| IHilton et al. 20T0l iTrarTet al. 20101) and 



steeper age gradients (IRosati et al. 20091 ) in clusters and over- 
dense regions at z ~ 1 - 1.6. 

The most successful method so far for finding high-redshift 
galaxy clus ters has been t hrough X -ray searches, archival or ded - 
icated (IRosati et al. 19981 iRomer et al. 200TllPierre et al. 20031) . 
Their depth is however constrained by the sensitivity of current 
observing facilities, which limits their effectiveness at higher 
redshifts. Colour selection techniques, on the other hand, 
can be used to efficiently search for clu sters up to z < 2 
(IGladders & Yee 20001 IWilson et al. 20081) and passive galaxy 
populations at even higher redshift (Kodama et al. 20071) . 
Finally, galaxy clusters might also be serendipitously discov- 
ered as spatial or redshift overdensities. The X-ray approach 
naturally selects massive and evolved structures, while red- 
sequence surveys are designed to search for a distinctive 
evolved galaxy population. Indeed, z > 1 X-ray and colour- 
selected structures are spatially compact and dominated 
by massive early-type galaxi es, as the i r local cou nterparts 
(IBlakeslee et al. 20031 IMullis et al. 20051 IStanford et al. 2 006, 
Pap ovich etal. 20101 ITanaka et al. 20101 |Henry et al. 20T0l 
Kur k et al. 2009b . In contrast, at z > 2 the search for clusters 
and their precursors has focused on finding overdensi- 
ties, often around radiogalaxies, of emission-line objects 
dFrancis et al. 19961 iPentericci et al. 19971 |Miley etal. 2 006, 
ISteidel et al. 20051 IQverzier et al. 2006b . in particular Lya 
emitters. Accordingly, high-density structures at z > 2 are char- 
acterised by a high level of star formation activity and mostly 
lack the extended X-ray emission and conspicuous early-type 
galaxy population typical of the evolved clusters. On the other 
hand, if the con centrati ons of Lya emitters reported at z > 3 
(ISteidel et al. 20001 IDaddi et al. 20091 IQverzier et al. 2008b are 
proto-cluster structures destined to evolve into the massive 
X-ray clusters observed at lower redshift, we should expect to 
find young yet already mature clusters at z ~ 2-2.5. However, 
whereas analogues to local massive clusters are known up to 
Z = 1.5-1.7, evolved early-type dominated and X-ray emitting 
clusters have not been found so far at earlier epochs, nor has 
an "intermediate" structure been observed at (or right after) the 
moment of quenching of star formation in its core galaxies. 



Here we present the discovery of CL J1449+0856, a con- 
spicuous galaxy overdensity at z = 2.07, dominated by mas- 
sive passively evolving galaxies and consistent with being the 
most distant X-ray detected galaxy cluster identified to date. In 
Section [2l we describe the target selection and photometric ob- 
servations. In Section (3j we present the high resolution imag- 
ing and morphological analysis of galaxies in the core and in 
Section |4l we discuss the spectroscopic observations and red- 
shift confirmation of the cluster. In Section we discuss the 
X-ray observations of the cluster and their analysis, while in 
Section [6| we compare CL J1449+0856 to other high-redshift 
structures. In Section [TJ we discuss its global properties and 
their implications for cosmology and Section [8] summarises our 
results. Unless specified otherwise, all magnitudes are reported 



in the AB system and we adopt a concordance cosmology with 
H = 70 km s" 1 Mpc" 1 , £l m = 0.3 and A = 0.7. 

2. Imaging and sample selection 

The structure was first identified, in archival Spitzer im- 
ages covering 342 arcmin 2 of the so-called "Daddi Field" 
(IDaddi et al. 2000b , as a remarkable overdensity of galaxies with 
IRAC colours [3.6] - [4.5] > at the position RA = 14h 49m 
14s and DEC = 8° 56' 21", indicating a massive structure at 
z > 1.5. Optical-NIR imaging data of this field was already avail- 
able from a multi-band survey in the B, R, I, z, and K s bands, the 
latter somewhat shallow. These data and their reduction are de- 
scribed in |Kong et al. 2006| and lDaddi et al. 20001 Between 2007 
and 2010, we obtained new deep imaging of this overdensity in 
the 7, /, H, and K s bands with MOIRCS on the Subaru telescope 
and in the / and K s bands with ISAAC on the VLT. For the 
purpose of studying galaxy morphology, we also obtained deep 
F160W imaging of the overdensity with NIC3 on Hubble and 
a shallow but high resolution K-band image using NIRC2 with 
adaptive optics on Keck. These two images were not used in the 
making of the photometric catalogue and their analysis is de- 
scribed in Section [3] In addition to the archival Spitzer/IRAC 
images, archival 24 [im data taken with Spitzer/MIPS were also 
available. We mention them for completeness but discuss them 
only briefly here, as they will be included in a future analysis 
of the galaxy population of the structure. Finally, X-ray observa- 
tions of the field were also available and are described in Section 
[5j Details of the imaging observations are given in Table [TJ 

The combined B, R, /, z, Y, J, H, and K s observations reach 
5cr limiting magnitudes of 26.95, 26.18, 26.03, 25.81, 25.64, 
25.47, 23.66 and 24.74 respectively. Catalogues were made 
for each band with SExtractor (Bertin & Arnouts 1996) using 
2"-diameter apertures and later merged. The final catalogue 
covers an area of 4' x 7, corresponding to the field of view 
of MOIRCS. As the galaxies with [3.6] - [4.5] > are better 
detected in the 7-band image than in the z and 5-band images, 
and since Y straddles the 4000 A break at z = 1.5, we do 
not rely on the traditional BzK criterion (IDaddi et al. 2004b to 
select for passively evolving galaxies at z > 1.5 but instead use 
Y - K s > 2, the expected colour of such a galaxy population. 
Out of 1291 objects in the combined catalogue with K s < 24.74, 
we find 114 red galaxies with Y — K s > 2. We note that only 1 1 
of those are detected at 5cr in B and 41 in z. 

The distribution of red galaxies shows a strong over- 
density at the same position as the IRAC-selected one. To 
characterise the overdensity, we created a number density map 
by dividing the MOIRCS field into a grid of sub-arcsecond 
cells and computing for each element the density estimator 
Z N = N/nr^, where r N is the distance to the Nth nearest galaxy 
with Y — K s > 2. We considered N = 3 - 7, which changed the 
angular resolution of the map but produced consistent results. 
The MOIRCS detector consists of two chips, one of which 
was centred on the overdensity and the other thus providing 
a low-density field. In the near-infrared images, the visible 
overdensity defines a 20"xl0" semi-axis elliptical area. We find 
that, in this 20" region centred on the structure, the mean density 
is 20cr above the field, with ~ 100 galaxies arcmin -2 in the 
overdensity versus five in the field. The density of this structure 
is thus similar to that report ed for the recen tly discovered 
galaxy cluster at z = 1.62 (Papo vTch et al. 2010| ). Fig. [T] shows 
two colour images of the overdensity, with contours represent- 
ing X-ray intensity and galaxy density, respectively, while Fig. 
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Fig. 1. RGB composite colour images of the 1.4'xl.4' field centred on the galaxy overdensity. The R channel of both images 
corresponds to the K s band; the G and B channels corresponds to the / and z bands in the left image and to the B and z bands in the 
right one. The B and z images were taken using the Suprime-Cam instrument on the Subaru telescope while the / and K s images are 
a composite of MOIRCS and ISAAC data, on the Subaru and VLT observatories. The BzK s image shows how the red galaxies are 
basically unseen at optical wavelengths. On the left, the white overlapping contours show the 1, 2, and 3<x significance levels of the 
diffuse X-ray emission in the XMM-Newton image, after subtraction of a point source seen in the Chandra image and smoothing 
with a 8" radius PSF (as described in Section 0. On the right, they show the 10, 20, and 30<x levels above the background galaxy 
number density, computed using the estimator. 
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Fig. 2. Surface density of galaxies with Y - K s > 2 in the whole MOIRCS field, in units of galaxies arcmin 2 as measured by the 
estimator. The black box delimits the 1.4'xl.4' field shown in Fig.[T] 
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Filter 


Central wavelength (p.m) 


Exposure (s) 


5cr limit (mag) 


Instrument 


Telescope 


Observation date 


B 


0.44 


1500 


26.95 


Suprime-Cam 


Subaru 


2003 Mar 5 


R 


0.65 


3600 


26.18 


Suprime-Cam 


WHT 


1998 May 19-21 


I 


0.80 


1800 


26.03 


Suprime-Cam 


Subaru 


2003 Mar 5 


z 


0.91 


2610 


25.81 


Suprime-Cam 


Subaru 


2003 Mar 4-5 


Y 


1.02 


17780 


25.64 


MOIRCS 


Subaru 


2009 Mar 15, 2010 Feb 7-8, 21 


J 


1.26 


9360 


25.47 


MOIRCS+ISAAC 


Subaru+VLT 


2007 Mar 10, Apr 5 


F160W 


1.60 


17920 




NIC3 


HST 


2008 May 11 


H 


1.65 


2380 


23.66 


MOIRCS 


Subaru 


2007 Apr 8 


K 


2.15 


1890 




NIRC2 


Keck 


2009 Apr 4 


K s 


2.20 


7800 


24.74 


MOIRCS+ISAAC 


Subaru+VLT 


2007 Mar 8, Apr 5 


IRAC 1 


3.6 


480 


23.85 


IRAC 


Spitzer 


2004 Jul 22 


IRAC2 


4.5 


480 


23.08 


IRAC 


Spitzer 


2004 Jul 22 


IRAC 3 


5.8 


480 


21.44 


IRAC 


Spitzer 


2004 Jul 22 


IRAC 4 


8.0 


480 


20.02 


IRAC 


Spitzer 


2004 Jul 22 


MIPS 24 


24 


480 


80 /Jy 


MIPS 


Spitzer 


2004 Aug 5 


0.5-10 keV 




80000 




EPIC-MOS 


XMM-Newton 


2001-2003 


0.5-8 keV 




80000 




ACIS 


Chandra 


2004 Jun 7-13 



Table 1. Details of the optical, near-IR photometric and X-ray observations. The B, R, I and z-band data have already bee n described 
in a previous paper (Kong et al. 2006). The XMM and Chandra data have also been described in other papers (IBrusa et al. 20051 
|Campisi et al. 2009] ). The raw IRAC and MIPS data were taken from the archive and reduced using MOPEX and custom scripts. 



|2]shows the number density of red galaxies in the MOIRCS field. 



3. High resolution imaging and galaxy morphology 

Very dusty star-forming galaxies at high redshift can have red 
colours similar to those expected for elliptical galaxies. While 
a complete morphological study of galaxies in the overdensity 
will be the subject of a future publication, we carried out here 
a first-order analysis to assess the nature of the red galaxies. 
We used if -band observations carried out with NICMOS-3 on 
Hubble, which provides the required sensitivity and spatial res- 
olution to unveil the morphology of the red galaxies, as well 
as AO-assisted ground-based imaging with NIRC2 on Keck. 
The NIC3 data were taken during seven orbits of Hubble. The 
individual frames were reduced using the NICRED pipeline 
(Magee, Bouwens & Illingworth 2007), which we found pro- 
vides a better calibration than the standard pipeline. In particular, 
the background noise varies less in the NICRED-reduced im- 
ages, which is critical to the morphological analysis as it allows 
us to recover the extended luminosity profile of faint galaxies. 
The individual frames w ere then combined using Multidrizzle 
(IKoekemoer et al. 20021) . The useful area of the NIC3 image, 
where the noise is low enough for the morphological analysis, is 
45"x50", covering the overdensity and its immediate surround- 
ings. 

The galaxies in this image were first matched with the 
catalogue and their surface brightness p rofiles modelled 
with a Sersic law (oc r l,n \ Sersic 1963) using GALFIT 
( |Peng et al. 2"002| ) and the single star in the NIC3 field as PSF. 
We fitted both the red (Y - K s > 2) and blue (Y - K s < 2) 
galaxies. For the latter, we used the BzK criterion to select for 
star-forming galaxies at z > 1.4: there are 352 sBzK- selected 
galaxies in the catalogue detected at 5cr in all three bands, of 
which ten are found in the NIC3 image. In that same field, 
we find 16 red galaxies. Fig. [3] shows the NIC3 image and 
the results of the Sersic fit: of those 16 red galaxies, five are 
unresolved (i.e. have an effective radius of one pixel or less; 
labelled "compact" in Fig. [3]) and the rest have n > 4, whereas 
three of the ten blue galaxies are unresolved, two have n > 3, 
and the rest n < 2. 

The NIRC2 image covers a similar area to the NIC3 image, 
but is rather shallow, having been taken during an unrelated ob- 
servation, and only the brightest galaxy cores are visible. It how- 



ever provides a good first-order verification of the Sersic mod- 
elling. Because of its combination of shallowness and very high 
resolution, the objects visible in the NIRC2 image either would 
be intrinsically bright or, in the case of faint galaxies like those 
in the overdensity, have a high surface brightness, similar to that 
of compact cores. Fig. |4] compares the NIRC2 and NIC3 im- 
ages. We find that, of the 13 red galaxies in the useful area of 
the NIRC2 image, eight are distinctly visible as very compact 
sources, supporting the results of the Sersic modelling. None of 
the sBzK galaxies rises above the noise of the NIRC2 image. 
The rest of the visible objects are obvious interlopers, low red- 
shift galaxies (colour- and angular size- wise) and a star. 

The galaxy overdensity is thus clearly dominated by 
spheroidal galaxies, as in local massive clusters. We note that 
more than half of these galaxies have r e > 2 kpc and thus appear 
less dense than previously studied passive galaxies at z > 2 
(IToft et al. 20071) . A thorough discussion of the stellar mass-size 
relation will be presented in a future paper. 



4. Spectroscopic observations and redshift 
determination 

We performed spectroscopic follow-up observations of the 
galaxies in and around the overdensity using FORS2 and 
VIMOS on the VLT and MOIRCS on Subaru. Blue, sBzK- 
selected galaxies were targeted around the structure's centre 
to within 10' (or 5 Mpc at z ~ 2; see Fig. 0. The VLT 
observations consisted of two FORS2 masks with 5 hours of 
integration each and one VIMOS mask with a 2.5h exposure. 
The 2D spectra of the individual runs w ere reduced using the 
standard pipeline ( Scodeggi o et al. 2005| ) and co-added. One- 
dimensional spectra were then extracted using the apextract 
tasks of the IRAF package. We estimated redshifts by first 
cross-correlating the observed spectra with a set of templates, 
including Ly man-break fS hapley et al. 2003 ), starburst, and 
star-forming (Kinn eyet al. 1996| ) galaxies. Using the rough 
redshift estimates given by the peaks of the cross-correlation 
function, we derived more precise redshifts from emission and 
absorption features (at z > 1.4, the FORS2 and VIMOS spectra 
cover the rest-frame UV) using the midlines task of the IRAF 
package. From the 41 FORS2 and 164 VIMOS slit spectra 
taken, we determined 109 secure redshifts. Their distribution 
shows a clear spike in the range z = 2 - 2.1, as shown in the top 
panel of Fig. [6] Assuming that this distribution is Gaussian, it 
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Fig. 3. Morphological properties of the galaxies in the field of the overdensity. Left, H band (F160W) NICMOS image of the cluster, 
in logarithmic greyscale. The image is a composite of frames taken with the NIC3 camera during 7 HST orbits, reduced using the 
NICRED pipeline and combined with Multidrizzle. Galaxies with Y — K s >2 are shown in red and galaxies with Y - K s < 2 in blue. 
Sersic indices of the best-fit model are indicated on top of each galaxy, except for those that are too point-like, which are labelled 
"compact". Right, morphologies of four representative galaxies from the H band image. For each galaxy we show, from left to right, 
the observed image, the best-fit Sersic model and the residuals image after subtraction of the model, all in the same logarithmic grey 
scale. 



Keck AO NIRC2 



0.16" - 1.3 kpc 



BCG? 




K-band,2.15 urn 



Fig. 4. Comparison of space and ground-based near-IR images of the field around the seemingly interacting galaxy triplet, taken 
with the NIC3 instrument on the HST (right) and NIRC2 on Keck using adaptive optics (left) and showing the compact cores of two 
of the three galaxies. In both images, the PSF is shown by a white circle. 



peaks at z = 2.07 and has a dispers ion of ~ 780 + 90 k m/s. Using 
the biweight estimator ( |Beers, Flynn & G ebhardt 1990}, we 
find z = 2.07 and 747 km/s respectively. These values of cluster 
velocity dispersion are comparabl e to that of star -forming 
galaxies in the nearby Virgo cluster (Bing geli et al. 1 987). With 



these values, we find 1 1 galaxies having spectroscopic redshifts 
within 2cr of the z = 2.07 peak. Some representative spectra of 
galaxies within the redshift spike are shown in Fig. [5] 



We also obtained spectra of several red galaxies in the 
near-IR with OHS/CISCO and MOIRCS on Subaru, but their 
faintness prevented us from measuring redshifts: while a 
stacked spectrum of the brightest member candidates shows 
some continuum, no absorption or emission features are seen, 
the latter down to typical limits of ~ 7 x 10" 17 erg s" 1 cm" 2 . At 
z ~ 2, this value corresponds to an unreddened star formation 
rate of < 10 M© yr" 1 and tends to support the conclusion that 
these objects are passively evolving stellar populations. 
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Fig. 5. Rest- frame UV spectra, taken with FORS2 and VIMOS on the VLT, of blue star-forming galaxies with redshifts within the 
peak of the distribution, centred at z = 2.07. The spectra shown here were rebinned with a bin size of two pixels. Prominent emission 
and absorption features are shown with dotted green lines. The dashed red line indicates the level of zero flux. For each object, the 
(uncalibrated) image from which the spectrum was extracted is shown at the top of the sub-plot. 



To complement the spectroscopic redshift information, 
which does not include the red galaxies, we estimated photo- 
metric redshifts from our BIRzY J HK S +IR AC photometry. Some 
galaxies in the catalogue have [5.8] - [8.0] > 0, suggesting that 
emission from an obscured AGN is contributing to the infrared 



SED and that the latter can therefore not be reproduced well by 
standard stellar population models. In particular, the galaxy in 
the overdensity with the highest IRAC excess is detected in the 
X-ray data, as discussed in Section[5j In these cases, we ignored 
the IRAC data for the purpose of deriving photometric redshifts. 
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Fig. 6. Redshift distributions of the galaxies in and surrounding the overdensity. Top: distribution of redshifts, determined from the 
emission-line spectra of blue (Y - K s < 2) galaxies taken with the FORS2 and VIMOS instruments on the Very Large Telescope. 
The insert shows a more detailed histogram for the range z = 2-2.15. For galaxies in the peak, the spectra sample the rest- frame 
~ 1200 - 2200A. Bottom: distribution of photometric redshifts of red (Y — K s > 2) galaxies within the central 20". These were 
estimated by comparing the BIRzYJHK s +IRAC spectral energy distribution to a set of spectral synthesis models computed from 
templates of half-solar, solar, and twice-solar metallicity and assuming two different types of star-formation history (exponentially 
declining and constant truncated). We did not include the effects of dust attenuation, as the colour, morphological, and spectral 
properties of the red galaxies point strongly towards them being dominated by passive stellar populations rather than dust-reddened 
younger stars. 



We compared the SEDs of the red galaxies to a range of model 
SEDs obtained from a set of Maraston 2005 stellar population 
synthesis templates, with three different star-formation histories 
(single burst, exponentially declining, and a constant star forma- 
tion rate with truncation) and metallicities ranging from half to 
twice the solar value. We did not include the effects of dust ex- 
tinction in the models, as the colour and morphology of the red 
galaxies, as well as the absence of emission lines in the spectra 
of those targeted with MOIRCS and OHS/CISCO, strongly sug- 
gest that they form a population of passively evolving systems. 
Fitting the SEDs with actively star-forming stellar populations 
and an arbitrary amount of dust results in a substantially worse 
X 2 - We also used the SED fit to estimate ages and stellar masses, 
assuming a bottom-light initial mass function (Chabrier 2003 ). 
Fig. |7] shows the SEDs and best-fit models of four red galaxies 
in the structure. We find that the distribution of the photomet- 
ric redshifts of red galaxies in the overdensity proper (i.e. within 
20" from the structure's centre; see Fig. [8]) narrowly peaks at 
z = 2.05 (Fig. [6]), with a scatter of ~ 0.07 that is fully com- 
patible with the expected accuracy for high-redshift ell ipticals 
(IDaddi et al. 20051 lllbert et al. 2"006l IMaraston et al. 20061) . We 
take the narrowness of both redshift distributions and their very 
similar peak values as confirmation that this structure, hereafter 
CL J1449+0856, is a real cluster or proto-cluster and we thus set 
its most likely redshift as z = 2.07. For the rest of the analysis, 
we consider as members of the structure galaxies with spectro- 
scopic or photometric redshifts within 2<x from the peaks of their 
respective distributions. We note that while z = 2.07 is the most 
likely redshift for this structure, having not yet measured any 



spectroscopic redshift for the red ellipticals and based on the 
widths of the spectroscopic and photometric redshift distribu- 
tions, we cannot exclude a slightly different redshift in the range 
2<z<2.1. 

Fig. [9] shows the spatial distribution of spectroscopic and 
photometric members compared to the size of the core and the 
field of view of MOIRCS. Table[2]gives the characteristics of the 
photometric members in the core of CL J1449+0856. 

4.1. Properties of the red galaxy population 

The colours of those red photometric members are significantly 
redder than those observ ed for established elliptical galaxies in 
z ~ 1.5 clusters (e.g. Papovich et al. 2010) and are consistent 
with a passive population at z ~ 2.1, as shown in Fig.[TOl We 
find a large colour scatter, however, indicating that the ellipticals 
have not yet settled into a tight red sequence characteristic of 
Z < 1.6 clusters. Indeed, their star-formation weighted ages (the 
mean age of the stars in the best-fit model) range from 0.6 to 2 
Gyr, with an average of 1.2 Gyr and average formation redshift 
of z ~ 3.5. This implies that, if truly passive, some of these 
galaxies are observed relatively shortly after the cessation of 
star formation and before the colour differences due to their 
different star-formation histories could be attenuated by the 
subsequent passive evolution. 

None of the elliptical members is much brighter than the 
others as to qualify for being a "brightest cluster galaxy" (BCG), 
the brightest of the red members being less than 0.5 mag more 
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Fig. 7. Spectral energy distributions of four representative red galaxies with photometric redshifts within 2<x of the peak of the 
distribution. The observed SEDs and errors are represented with blue open circles. Upper limits at l<x are shown by arrows and 
the best fit template is shown in orange. The corresponding redshift and formation redshift, at which half of the stellar mass was in 
place, are given in the upper right corner of each sub-plot. The integrated template fluxes in the Suprime-Cam, MOIRCS and IRAC 
bands are shown by red squares. 



luminous than the second and third brightest members. However, 
as seen in Figs. [3] and |H a group of three galaxies near the centre 
of the cluster appear to be very close to each other. If they all are 
at the same redshift of z = 2.07, they are separated by 5.5 to 13 
kpc in projection and thus likely interacting. As noted in Section 
[3] Keck adaptive optics observations with NIRC2 reveal in two 
of them a very compa ct (~ 1 kpc) co re typical of high-redshift 
elliptical galaxies dDaddi et al. 2005b . The combined flux of 
the three components is ^ab, rest-frame ~ 20, consistent with 
the /T-band-redshift relation for BCGs ( |Whiley et al. 2008) , 
when extrapolated to z = 2.07. This suggests that we might be 
witnessing the early stages of th e assembly of a bri ghtest cluster 
galaxy through merging (De Lucia & Blai zot 20071) . 



We note however that some emission is seen at 24 /mi in and 
around the cluster centre. Several sources are detected at 5cr in 
the core (Fig. [TT]), some of them associated with red galaxies, 
including the AGN and the "proto-BCG" triplet. The PSF-fitted 
fluxes of these objects are 100-150 /Jy which, if due to star 
formation, would correspond to ULIRG-like luminosities of 
L IR ~ 1.5-2xl0 12 L© and imply very active starbursts with 
star formation rates of > 150 M© yr" 1 , in apparent contradiction 
with the elliptical-like morphology and compact cores. In the 
case of the "proto-BCG" galaxy group, the MIPS PSF is too 
large (6") to determine which of the three galaxies is a 24 jim 
emitter. On the other hand, the 24 fim emission (which at z ~ 2 
corresponds to 8 yum rest-frame) might be due to extremely 



obscured AGN activity (Fio re et al. 20091 IDaddi et al. 2 007b K 
an interpretation supported by the presence among the MIPS 
24 jim sources of the Chandra-detected AGN. Stacking the 
Chandra data at the position of the five other 24 [im detected 
red galaxies within the central 1', using a 5" aperture, we find 
an excess of photons in the hard band at 2.3<x significance with 
respect to the distribution of counts in a thousand samples of 
five random background (i.e. chosen to be at least 5" away 
from the nearest K s selected object) positions each. As reported 
in Section we find no such excess in the soft band. As we 
might include X-ray emitting but ^-undetected sources in the 
background positions, the significance level is likely underes- 
timated. However, distinguishing the relative contributions of 
AGN activity and star formation to the mid-IR flux is not within 
the scope of this work and would require far-IR data. 



5. XMM-Newton and Chandra imaging 

The unambiguous signature of an evolved cluster is the 
X-ray emission from the ICM, as it implies a deep and 
established potential well. To check for emission from a 
diffuse atmosphere, we looked for extended emission in 
deep X-ray observations available in the field with both the 
XMM-Newton and Chandra telescop es totalling 80 ks each 
(IBrusa et al. 20051 [Campisi et al~2 009). A detection was found 
in both soft-band (0.5-2 keV) images: the Chandra observation 



8 



Gobat et al.: An old cluster in the young Universe 



scole ot z = 2.07 (Mpc) 
-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 




222.320 222.315 222.310 222.305 222.300 

RA (deg) 



Fig. 8. RGB composite image of the 1 .4' x 1 .4' field centred on the CL J1449+0856, as in Fig.Q] The R, G, and B channels correspond 
to the K s , J, and Y bands, respectively. The Y, J, and K s photometry straddles wavelengths close to the 4000 A break at z ~ 2 and 
is useful to appreciate likely age variations among potential cluster members. The red arrows show the positions of red galaxies 
with photometric redshifts within 2<x of the peak of the distribution and the blue ones the position of galaxies with spectroscopic 
redshifts within 2<x of the peak. The green square gives the position of the AGN seen with Chandra. The white ellipse shows the 20" 
semi-major axis region that corresponds roughly to the overdensity and from which were selected the red galaxies shown in Figs. [7] 
andHOland in Tabled 
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Table 2. Characteristics of the red galaxies in the overdensity: coordinates, photometric redshift, Y - K s colour, ^-band total 
magnitude, and distance from the cluster centre in arc minutes and kiloparsecs. The first ten galaxies are the red members, with 
1.9 < Zphot < 2.2. The bottom six objects, which are not formally inside the photometric redshift peak, might still belong to the 
structure: e.g. galaxy 51 13, which is the Chandra-detected AGN, a red spheroidal (nsersic = 6) whose optical-NIR SED likely suffers 
from AGN contamination. 



reveals a ~1" point source at the position of one of the red 
galaxies, while the soft X-ray emission seen by XMM is more 
extended than the 6" instrument PSF. To assess the significance 
of this extended emission, we fitted the XMM emission with a 
PSF profile at the position of the Chandra source and analysed 
the residual image. After subtraction of the point source, we 
found a residual X-ray emission in the soft band at the 3.5<x 
level on scales of 20-30", three to five times more extended than 
the PSF of XMM-Newton. The excess flux over the background 



is 47 + 13 photons, the error including systematic uncertainties 
due to the point source subtraction. The total flux of this 
extended emission in the range 0.5-2 keV and over a 16" radius 
is ~ 9.3 x 10" 16 erg s" 1 , consistent with the presence of hot 
ICM, typical of a "relaxed" cluster. This is thus faint emission, 
yet detected at a higher signific ance level than what rep orted for 
the z = 1.62 cluster (ITanaka et al. 20101 [Papovich et al. 20T0| . 
As a consistency check, we carried out the same analysis on 
the Chandra data. While we do not detect any residual extended 
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Fig. 9. Spatial distribution of blue spectroscopic members, centred on the cluster. Galaxies with spectroscopic redshifts within 2cr 
of the z = 2.07 peak are shown by blue open circles. The positions of galaxies with z p hot within 2<x of the z = 2.05 peak are shown 
by red crosses. The orange rectangle shows the field of MOIRCS and the 1, 2, and 3<x significance levels of the extended X-ray 
emission (as shown on Fig. [T] and discussed in Section© are drawn at the centre of the plot. 



emission in the Chandra image, the implied upper limit is 
nevertheless consistent with the XMM result. Fig. [12] shows the 
XMM and Chandra detections and the extended emission after 
subtraction of the Chandra point source. 

We emphasise that we kept the X-ray analysis simple and 
robust. We di d not apply so phisticated wavelet-based detection 
techniques (Finogu eliov et al. 2006] ), but only a standard back- 
ground subtraction, point-source subtraction, and photometry 
over a 16" radius. For visualisation purposes only, we smoothed 
the data with a 8" radius PSF to determine the X-ray contour 
levels shown in Fig. [T] For an accurate subtraction of point 
source X-ray emitters, we used the deep Chandra observations 
of the field, exploiting the fact that Chandra has a much better 
spatial resolution than XMM. On the other hand, we note that, as 
the Chandra data were taken one to three years after the XMM 
observation s (Table [p, AGN var iability could have affected 
our result (Papada kis et al. 2008] ). Our measurement should 
however be robust against variability for two strong reasons: 
first, the fact that the XMM residual emission is recovered 
over an extended area, after point source subtraction, implies 
that it is not an effect of variability. Furthermore, with the 
position of the AGN known, we fitted its flux using the XMM 
data itself only and found results fully consistent with what 
is observed in Chandra, excluding substantial variability. We 
also considered the possibility that the extended X-ray emission 
might actually arise from several faint AGNs, not individually 
detected by either XMM or Chandra and spatially dispersed in 
the overdensity. While detailed study of the AGN content in this 
cluster is deferred to a future paper, we carefully investigated 
which of the galaxies in the structure might be hosting an AGN 
on the basis of either a hard X-ray (2-8 keV) detection, the 
presence of a rising IRAC SED, or of mid-IR excess emission 
(Daddi et al. 2007b). For the latter, we used the archival Spitzer 



MIPS 24/zm imaging of the field, which reaches a 5<x detection 
limit of about 80 //Jy. We find that there could be up to 4-5 
AGNs in the cluster (as shown in Fig. [TT]), in addition to the 
Chandra point source. However, no soft X-ray photons were 
detected by Chandra at their position. Finally, we considered 
the possibility that the X-ray emission might come from a lower 
mass foreground (e.g. z < 1) galaxy group. Such a structure 
would be more loose and not immediately apparent against the 
backdrop of the main galaxy overdensity but should appear 
as a distinct peak in redshift space. To check for foreground 
structures, we estimated photometric redshifts for Y — K s < 2 
galaxies without spectroscopic redshifts. For these objects, 
which include star-forming galaxies as well as low redshift 
ellipticals, we used Coleman, Wu & Weedman 1980 templates, 
with the addition of a 100 Myr constantly star-forming model 
computed from Marast on 20051 templates, and included dust 
extinction up to E(B - V) = 1 . We find no significant secondary 
structure in the redshift distribution of the 33 galaxies (3 
spectroscopic and 30 photometric redshifts) within the lcr 
confidence region of the extended X-ray emission. Three blue 
galaxies with z P hot ~1.7 might be associated but their projected 
centre of mass is offset by 15" with respect to the X-ray centroid. 



We conclude that, all in all, the only likely explanation is 
that the extended X-ray emission seen by XMM is due to an 
ICM present in the structure's potential well. The presence of 
an X-ray atmosphere and an evolved galaxy population is con- 
sistent with CL J1449+0856 being not a forming proto-cluster 
but an already mature galaxy cluster comparable to the massive 
structures observed at z < 1.6. 
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Fig. 10. Colour-magnitude diagram of the objects within 20" of the cluster centre. Blue galaxies are shown by blue crosses and 
red galaxies by red squares. Red galaxies with a photometric redshift within 2<x of the peak of the distribution are indicated by 
filled squares, the bars showing t he error in the colour. The greyscale shaded map shows the distribution of #z^-selected galaxies 
in the GOODS-South field (Dad di et al. 2007 at . Colour limits at 3<x are shown by arrows. The typical error in the Y - K s colour 
and K Sttot magnitude is shown in the bottom left. The solid and dashed lines show the expected colour and magnitude of stellar 
population synthesis models (Kodama & Arimoto 1997), assuming that star formation begins at z = 10 and z = 3.5, respectively. 
The composite flux and colour of the "proto-BCG" galaxy assemblage are marked by a green star. 



> 2 have 



Miley et al. 2006 



6. CL J1 449+0856 and z ~ 2 structures 

While several overdensities of Lya emitters at z 
been reported (Penteric ci et al. 19971 
ISteidel et al. 2"005l IQverzier et aL2006l CL J1449+0856 
is different in some key aspects because of the evidence of 
extended X-ray emission from an intra-cluster medium and a 
centre occupied by old passively evolving early-type galaxies. 
These features make CL J1449+0856 much more mature than 
the high-redshift proto-clusters and a unique case among the 
z > 2 structures. 



6.1. MRC 1138-262 



MRC 1138-262 has been observe d with Chandra and 
a number of AGNs were detected (Pente ricci et al. 2002b 
as well as diffuse emission centred on the radio galaxy 
dCarilli et al. 2002b . The latter is however clearly associated 
with the central AGN, which contributes to more than 80% 
of the total flux, aligned with the radio lobes, suggesting 
that it is due to inve r se Compton scattered CMB photons 
(ICelotti & Fabian 20041 |Finoguenov et al. 2010| ), and overall 
quite different from the typical extended emission from an 
intra-cluster medium ( Carill i et al. 20021) . On the other hand, 
the only AGN in the X-ray emission of CL J1449+0856 does 
not contribute more than -50% of the total observed flux. The 
extended emission itself is not centred on any particular galaxy, 
but clearly spatially associated with the galaxy overdensity, 
fully consistent with a young cluster atmosphere. 



Among the large galaxy overdensities at z > 2, 
that around the massive radio galaxy PKS 1138-262 
(IPentericci et al. 19971 |Miley et al. 2006] ) at z = 2.16 is the 
one that can be most readily compared to CL J1449+0856. 
Located at a similar redshift, and thus at the same epoch in the 
history of the Universe, it is a massive structure characterised 
by a giant radio galaxy and a host of significantly less massive 
star-forming "satellite" galaxies (Mile y et al. 2006| ). It has been 
extensively studied photometrically as well as spectroscopically 
and is understood to be a structure still in its formation phase, 
i.e. a "proto-cluster". We emphasise that there are significant 
differences indicating that CL J1449+0856 is a very different 
type of structure, observed at a substantially more advanced 
evolutionary stage than MRC 1138-262. 



The proto-cluster MR C 1138-262 was ident ified as an over- 
density of Lya emitters (Penteri cci et al. 200 0) and the radio- 
galaxy itself is emb edded in a giant emission-line nebula 
( Pente ricci et al. 19971) . The core of the proto-cluster is entirely 
dominated by star-forming galaxies, with a few redder galaxies 
located at the outskirts of the Lya halo dHatch et al. 20091). O f 
the latter, two have detected Ha emission ( Doh erty et al. 2010| ), 
indicating that they are actually actively star-forming galaxies. 
Red galaxies have been found in a broader field surrounding the 
radio-galaxy and Lya overdensity, a subset of which have SEDs 
and morphologies consistent with a passive population, but the 
claimed overden sity is not concentrated in a putative cluster core 
(IZirm et al. 2008b , mostly avoiding the central region around the 
radio-galaxy. In contrast, the structure of CL J1449+0856 is very 
different. Simply looking at colour images shows the striking 
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Fig. 11. ^-band image of the 1.4'xl.4' field centred on the galaxy cluster, as in Figs. \T\ and [8j the green contours showing the fluxes 
of the sources detected at 5<x or more in the MIPS 24 jim image. The red arrows show the positions of red (Y - K s > 2) galaxies 
with rising IRAC fluxes ([5.8] - [8.0] > 0), suggesting the possible presence of obscured AGN activity. 




Fig. 12. Left and middle: soft (0.5-2 keV) X-ray images of the field taken with XMM-Newton and Chandra (80ks each), showing, 
respectively the diffuse emission and a point source, an active galactic nucleus corresponding to one of the red galaxies. Right: 
signal-to-noise map of the XMM residual image, after subtraction of the point source and smoothing with a FWHM of 16". 

difference between our clust er and MRC 1138 -262 (compare 
Fig. 1 in this paper to Fig. 1 in Hatch et al. 2009, which displays 
a comparable field of view). Specifically, no analogue to the Lya 
halo of MRC 1138-262 has been observed in CL J1449+0856. 
None of the red galaxies in the core is seen in the rest- frame UV 
and the space between the galaxies is void of light, save for the 
diffuse X-ray emission. And whereas the core of MRC 1 138-262 
is dominated by star-forming galaxies, the core population of CL 
J1449+0856 is constituted of red, morphologically early-type 
galaxies. This supports our conclusion that, if CL J1449+0856 
experienced such a star-forming phase as the proto-cluster is un- 
dergoing, it lies in its past and that the red galaxies are now pas- 
sively evolving. 



Finally, MRC 1138-262, having been discovered as a pow- 
erful FRII source ( Pente ricci et al. 1997b , is in practice selected 



over the whole sky and hence it is very difficult to assess its cos- 
mological relevance. 

6.2. JKCS 041 

lAndreon et al. 2009l claimed the discovery of an evolved, colour- 
selected galaxy cluster at z ~ 1.9, with extended X-ray emis- 
sion detected by Chandra. From the X-ray luminosity, the au- 
thors derive a total mass well above 10 14 M©. The redshift of 
the structure, JKCS 041, was estimated from photometric red- 
shifts only: red-sequence galaxies were observed with FORS2 
but proved too faint and no spectroscopic information was thus 
provided to support this identification. The redshift identification 
of this cluster has since been disputed usi ng spectroscopic red- 
shifts from the VIMOS-VLT Deep Survey (ILe Fevre et al. 2005b 
and new multi-band photometric redshifts; these data show that 
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JKCS 041 is not a single high-redshift cluster but a superposition 
along the line of sig ht of at least two r ich galaxy structures, at 
z ~ 1.1 and z ~ 1.5 ( jBielby et al. 2010) . While a z ~ 2 structure 
could still lurk in the background, this result suggests that the 
X-ray emission originates at z < 1.5. 

7. Structure and mass of CL J1 449+0856 

With an estimated redshift of z = 2.07 for the cluster, we derive 
an X-ray luminosity of L x (0. 1 - 2.4 keV) = (7 + 2) X 10 43 erg s" 1 . 
On the basis of established L x - M correlations 
(Leauthaud et al. 2QT0b . the luminosity corresponds to a 
total mass of M 2 oo = (5.3 ± 1) x 10 13 M Q , comparable 
to that of Virgo. The corresponding virial radius would be 
^200 ~ 0-37 Mpc, consistent with the scale of the observed 
XMM emission (~ 0.4 Mpc, as shown in Fig. []]), and a 
temperature of kT ~ 2 keV, below the detection limit of 
current SZE facilities ( Vanderlinde et al. 2010b . We derive 
a consistent estimate for the total mass from the total stel- 
lar mass of the red member galaxies in the central 20", 
M* t = 4.9 x 10 11 M©. Using t he locally calibrate d stellar 
mass-to-halo mass conversion of IMoster et al. 20101 we find 
a halo mass of Mh a i = 8.0 x 10 13 M Q and a corresponding 
virial radius of 0.42 Mpc (assuming Mh a i =M2oo)- We note that 
both estimates of the virial radius are comparable to the size of 
the overdensity from which we selected the red galaxies, 20" 
corresponding to ~ 0.17 Mpc proper at z = 2.07. Furthermore, 
assuming that the core is in the first order spherical, at z = 2.07 
a sphere of radius 0.17 Mpc with a mass of 8 x 10 13 M© has 
a mean density of ~ 3000p c (where the critical density at 
z = 2.07 is p c = 8.63 x 10" 29 g cm -3 ), significantly above 
the usual threshold of ~ 178p c for virialisation. Although 
this estimate has substantial uncertainties, this might suggest 
that the core of CL J1449+0856 has already collapsed and (at 
least partially) virialised. Both total mass estimates might be 
regarded as lower limits, as the L x - M relation assumes that 
the intra-cluster gas has fully reached its virial temperature, a 
state that CL J1449+0856 may still be approaching, and our 
census of massive galaxies is incomplete as we considered only 
the red members in a small central region, without including 
star-forming members or galaxies outside the immediate core. 
The agreement of the two estimates, however, suggests that 
the faint X-ray emission is reliable. It would seem to imply 
that not only the M* - M relation holds to z ~ 2 and but that 
this cluster is, to the first order, already on the L x -M correlation. 

With these halo mass estimates, we can now use the 
results of numerical simulations from the literature to attempt 
a prediction of the future growth of CL J1449+0856 and its 
final mass at z = 0. Based on the mass assembly histories 
of halos in the Millennium and Millennium-II simulations 
(IFakhouri et al. 2010b . we estimate that CL J1449+0856 would 
reach a mass of 0.9 - 1 .5 x 10 14 M at z = 1 .5 and 4.9 - 8.2 x 10 14 
at z = 0. These values are well within the range of cluster masses 
at these respective redshifts. At z = 1 - 1.5, CL J1449+0856 
would be about one fourth as massive as the most ma ssive 
z > 1 clusters (IBrodwin et al. 20101 IRosati et al. 2009b and 
reach a z = m ass comparable to that of the Coma cluster 
(IKubo et al. 2007b . 

Assuming Gaussian in itial conditions and concordance cos- 
mology with erg = 0.8 dVan derlinde et al . 2010b . the probability 
of finding a dark matter halo with z > 2.07 and a mass greater 
than 5 x 10 13 M© in the survey area of 400 arcmin 2 is 3.3 x 10" 2 . 



Allowing some time before the dark matter halo formation and 
observation (Jimenez & Verde 2009), and the gas to settle in the 
deep potential well, would further lower this probability, e.g. 
down to 1.0 x 10" 2 for Zf = 2.30 (0.3 Gyr before the observa- 
tions). Lowering erg to 0.77 would further reduce these numbers 
by roughly a factor of 3. Although based on a single object de- 
tection, our finding suggests an excess of massive high-redshift 
structures, consistent with the constraints derived from the ex- 
istence of a very mass ive cluster at z = 1.39 dJee et al. 2 009, 
IJimenez & Ve rde 2009). However, we note that no obvious sim- 
ilarly high-redshift candidate structure is detected in the 2 deg 2 
COSMOS field ( Tanaka et al., in prep^} Salvato et al., in prep!] ). 



8. Summary 

We have discovered a remarkable structure whose properties 
are consistent with it being a mature cluster at z = 2.07. This 
structure was selected as an overdensity of sources with IRAC 
colours satisfying [3.6] - [4.5] > 0. Deep follow-up observa- 
tions with Subaru and the VLT revealed a strong overdensity of 
galaxies with Y-K s colours consistent with a passive population 
at z > 2. We also obtained high resolution HST/NICMOS and 
Keck AO images, which revealed that the red galaxies have 
elliptical-like morphologies and compact cores. 
From the VLT VIMOS and FORS2 spectra of sBzK galaxies 
around the core, we estimated a redshift of z = 2.07 for the 
structure. We estimated photometric redshifts from the 12-band 
SEDs of the red galaxies, whose distribution peaks at z = 2.05, 
in agreement with the spectroscopic redshifts. 
Using XMM-Newton and Chandra observations of the field, we 
found an extended soft X-ray emission at the 3.5<x confidence 
level, at the position of the galaxy overdensity. The observed 
X-ray luminosity and the galaxy mass content of the core imply 
a total halo mass of 5-8xl0 13 M . 

Our results show that virialised clusters with detectable X- 
ray emission and a fully established early-type galaxy content 
were already in place at z > 2, when the Universe was only 
~ 3 Gyr old. While it took us several years of observations to 
confirm this structure, upcoming facilities like JWST and future 
X-ray observatories should be able of routinely find and study 
similar clusters, unveiling their thermodynamic and kinematic 
structure in detail. The census of z > 2 structures similar to CL 
J1449+0856 will subject the assumed Gaussianity of the primor- 
dial density field to a critical check. 
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